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ABSTRACT 

In conventional models of galactic and cluster cooling flows widespread cooling (mass dropout) is 
assumed to avoid accumulation of unacceptably large central masses. This assumption is often criticized 
because the observed mass of cooled gas and newly formed stars are much less than the cooling rates 
imply. Moreover, recent XMM observations have failed to find spectral evidence for locally cooling gas. 
This remarkable discovery has revived the notion that cooling flows do not cool in the manner previously 
assumed, but are instead heated by some process such as an intermittent or low-level AGN in the central 
galaxy. To explore this hypothesis, we consider the gas-dynamical consequences of heated galactic cooling 
flows under many different heating scenarios without specifying the detailed physics of the heating 
process. We have been unable to find a single acceptable heated flow in reasonable agreement with well- 
observed hot gas temperature and density profiles. Even for poor fits to the observations, these models 
require finely tuned heating scenarios. Idealized flows in which radiative cooling is perfectly balanced by 
global heating are grossly incompatible with observations. Flows heated by central feedback often result 
in eposidic cooling that is associated with quasi-cyclic changes in the hot gas density profile; at present 
there is no observational evidence for this. Paradoxically, when cooling flows are centrally heated (or 
pressurized), they experience spontaneous non-linear compressions that result in spatially widespread 
cooling instabilities. Therefore, spectral evidence for cooling gas is expected in either traditional cooling 
flows with ad hoc mass dropout or in heated cooling flows without the dropout assumption. 

Subject headings: galaxies: elliptical and lenticular - galaxies: active - galaxies: cooling flows - 
galaxies: radio - X-rays: galaxies - X-rays: galaxy clusters 



1. INTRODUCTION 

Early results from the XMM and Chandra X-ray tele- 
scopes have confounded and intensified long-simmering 
controversies regarding cooling flows in elliptical galaxies 
and galaxy clusters. In the absence of detailed spectro- 
scopic data, it had been assumed that the virialized hot 
gas in these systems cools over a large volume as the gas 
slowly flows inward, depositing several to many 100 A'Iq 
yr^^. Since optical evidence for normal star formation is 
limited, it is usually assumed that the cooled gas forms 
preferentially into stars of low mass. Indeed, there is some 
evidence for cooled gas in Ha -I- [Nil] emission in cool- 
ing flows (e.g. Heckman et al. 1989; Caon, Macchetto & 
Pastoriza 2000) and for populations of younger stars (Mc- 
Namara 1997; Cardiel et al. 1998). Gas dynamical mod- 
els with spatially distributed mass dropout (cooling) have 
been devised that agree with most observations. Many 
X-ray observers have claimed that cooling flows can be 
understood if mass drops out so that the remaining mass 
flow varies as M{r) cx r (e.g. Thomas, Fabian, & Nulsen 
1987). 

The unexpected result from the Reflection Grating 
Spectrometer (RGS) on XMM is that the multiphase spec- 
tral signature of gas cooling to low temperatures has not 
been seen in Abell 1835 (Peterson et al. 2001) or Abell 
1795 (Tamura et al. 2000). XMM/EPIC spectra of 



M87/Virgo also fail to reveal evidence of cooling at the 
expected rate (Bohringer et al. 2001a, b, c; Molendi 
& Pizzolato 2001; Matsushita et al. 2002), in contra- 
diction to the long-standing observation with Einstein 
FPCS by Canizares et al. (1982), recently confirmed by 
Buote (2001). High resolution XMM/RGS spectra of el- 
liptical galaxy NGC 4636 show strong emission from the 
OVIII Lya line at 0.654 keV characteristic of hot gas at 
T > 6 X 10^ K, but the OVII resonance line at 0.574 keV, 
expected from cooling gas at T 2 x 10^ K, is strangely 
absent (Xu et al. 2001). On the other hand, FUSE ob- 
servations of NGC 4636 detect OVI emission from gas at 
T ~ 3 X 10^ K indicating a mass cooling rate M near (but 
less than) that predicted by conventional multiphase cool- 
ing flow models (e.g. Bregman et al. 2001). Nevertheless, 
as a result of the XMM data, an increasing number of the- 
orists and observers are doubting the reality of traditional 
cooling flows and radically modified alternative theoretical 
models have been proposed (Fabian et al. 2001; Morris & 
Fabian 2001). 

One of the most consistent and (initially) appealing sug- 
gestions is that the gas in these so-called cooling flows is 
not cooling by radiation losses because it is continuously 
or intermittently heated by energy released near supermas- 
sive black holes in the cores of central E or cD galaxies. 
This seems reasonable since almost every giant E or cD 
galaxy is known to be a source of AGN activity at some 
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level, ranging from the barely observable to the most pow- 
erful known radio galaxies. It has been suggested that the 
hot thermal plasma could be heated by AGN energy in 
the form of shock waves, Compton scattering or by dis- 
sipation of the energy in relativistic particles. Certainly 
some energy transfer to the gas must be occurring, and 
the AGN is likely to operate in a feedback mode, supply- 
ing energy whenever gas cools at the center of the cooling 
flow and fuels the black hole. Enough AGN energy is avail- 
able to significantly influence the flow, provided it can be 
judiciously transfered to the gas. 

Although the heating hypothesis is appealing, we show 
in the discussion below that unforeseen complications rule 
out many - and perhaps all - heating scenarios that af- 
fect the gas flow signific;antly. Our objective here is to 
study the time dependent evolution of cooling flows that 
are heated in a variety of different ways. We do not con- 
sider the detailed and uncertain microphysics of the heat- 
ing process. Instead, we assume that heating occurs ac- 
cording to various space and time dependent recipes as 
described by a few adjustable parameters. For success, a 
heated flow must agree reasonably well with the density 
and temperature distributions observed in cooling flows. 
In addition, this agreement must be achieved with a wide 
range of heating parameters so that it is not necessary to 
"fine-tune" the heating process, which may be inherently 
stochastic. Furthermore, the mass of gas accumulated at 
or near the center of the flow must be consistent with the 
masses of central black holes and stellar mass to light de- 
terminations from central velocity dispersions. Finally, to 
be fully compatible with the lack of evidence for intrin- 
sic multiphase gas in XMM X-ray spectra, heated cooling 
flows shoiild have little or no distributed cooling. 

We focus here on galactic or galaxy group scale cool- 
ing flows associated with a single central massive elliptical 
galaxy. Indeed, all giant E galaxies are thought to have 
formed by mergers of smaller galaxies in group environ- 
ments. As a result, most luminous E galaxies are sur- 
rounded by extended X-ray halos confined by the group 
potential, containing more hot gas than can be accounted 
for by stellar mass loss in the central galaxy. Within the 
central galaxy, however, new gas is continuously supplied 
by stellar mass loss so the absence of observed cooling flow 
dropout is consequently more of a problem. Since larger 
cluster cooling flows invariably flow into a central E or 
cD galaxy, our arguments generally apply for these larger 
flows as well, although more central heating would be re- 
quired to decelerate the inflow velocity in cooling in cluster 
flows. 

We deliberately choose the simplest possible heating sce- 
narios; if these do not give satisfactory results, it is unlikely 
that adding details will help. We warn readers that most 
or all of the heated cooling flow models described here are 
found to be unacceptable by the criteria that we have set. 

It is particularly interesting that multiphase cooling over 
a large region of the flow is a frequent and natural conse- 
quence of heating or additional nonthermal pressure near 
the center of the galaxy. Centrally heated or pressurized 
gas decelerates (or reverses) the normal cooling inflow in 
the galactic core. At the radius where normally inward 
flowing gas from large radii encounters stagnant or out- 
flowing gas in the core, a nonlinear compression develops. 



resulting in thermal instabilities. Since this cooling in- 
stability is inherently nonlinear, it was not found in linear 
studies (e.g. Balbus & Soker 1989). However, the expected 
multiphase spectral signature for such cooling gas is not 
yet reported by XMM observers. 

AGN heating in cooling flows is not a new idea. Before 
we describe our own models, we review in the following 
section some of the heated models that have already been 
proposed. 

2. A BRIEF HISTORY OF HEATED COOLING FLOWS 

2.1. Possible Heating Mechanisms 

2.1.1. Thermal Conduction 

The positive temperature gradient {dT/dr > 0) in the 
central regions of cooling flows is due to radiative cooling 
and, on galactic scales, to mixing with cooler stellar mass 
loss. Since the heat flux opposes the thermal gradient, 
it has been suggested that heat could be conducted in- 
ward toward the center of the flow, balancing the radiative 
cooling there and arresting or retarding the global evolu- 
tion of the cooling flow and its mass dropout (Stewart, 
Canizares, Fabian & Nulsen 1984; Bertschinger & Meiksin 
1986; Meiksin 1988; Bregman & David 1988; Loewenstein, 
Zwcibcl & Begelman 1991; Brighenti & Mathews 2002). 
Unfortunately, these attempts to avoid radiative cooling 
have been generally unsuccessful because conductive heat 
flux and radiation losses are important in different flow re- 
gions. The Spitzer conductivity depends strongly on tem- 
perature {k oc T^/^) and is often very small at the centers 
of cooling flows where radiative losses (oc n^) are greatest. 
Furthermore, when the temperature is low everywhere in 
the flow, T < 1 keV, as in galactic scale cooling flows, ther- 
mal conductivity is small at every radius. Thermal con- 
duction can be reduced by tangled magnetic fields, which 
may explain the existence of "cold fronts" in some cluster 
cooling fiows (Ettori & Fabian 2000, but also see Rosner & 
Tucker 1989; Narayan & Medvedev 2002). For galaxy clus- 
ters, the inward conduction of thermal energy can slow or 
stop the cooling flow, but the conductivity must be care- 
fully adjusted to fit a particular observed flow. Generally 
the results are very sensitive to the conduction reduction 
factor and are cluster-dependent, a form of fine-tuning that 
invalidates most models. 

2.1.2. Heating by Cosmic Rays 

Radiofrequency emission from cosmic ray electrons is 
observed in the galactic centers of most or all cooling flows. 
If the less visible proton cosmic ray component is scaled 
as in the Galaxy, Coloumb collisions could produce sub- 
stantial heating (Rephaeli 1987; Rephaeli & Silk 1995). 
Hints that this type of heating could be important are 
just now emerging. By combining the most recent X-ray 
and radio data from the X-ray holes in NGC 1275 (Perseus 
cluster), Fabian et al. (2002) find evidence for an energy 
density of suprathermal particles that exceeds that of the 
fast electrons by (180 — 500)/, where / is the (unknown) 
volume filling factor of the electron component. In addi- 
tion, the flattened gas cores at r < 700 pc in bright Virgo 
E galaxies strongly suggests the presence of an additional 
non-thermal pressure component which could be due to 
(non-radiating) cosmic rays or magnetic flelds (Brighenti 
& Mathews 1997a). More commonly, cosmic ray heating is 
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assumed to be due dissipation of Alfvcn waves excited by 
relativistic electrons or protons streaming along field lines. 
This possibility has been widely discussed in previous work 
with somewhat inconclusive results (Lerche 1967; Kulsrud 
& Pearcc 1969; Wentzcl 1969; Skilling 1971; Scott et al. 
1980; Tucker & Rosner 1983; Rose et al. 1984; Rosner & 
Tucker 1989; Begelman & Zweibel 1994). 

2.1.3. Inverse Compton Heating by hard AGN Radiation 

Compton heating by isotropic quasar radiation (Bin- 
ney & Tabor 1995; Ciotti & Ostriker 2001) is effective if 
the Compton temperature Tq ~ {hi/) /4 exceeds the virial 
temperature, typically Tyir ^ 1 keV 10^ K) for galac- 
tic cooling flows. Gamma ray observations with EGRET 
since 1991 have discovered ^ 50 flat spectrum blazars (BL 
Lac AGN and OVV type QSOs) in which Tc ~ 10« - 10^ 
K due to the dominance of gamma rays in the energy out- 
put (von Montigny et al. 1995;. Impey 1996; Hartman et 
al. 1999). If all normal elliptical galaxies emitted such ra- 
diation isotropically following short-lived accretion events, 
cooling flow gas near the AGN, when heated to Tc ^ Tyi,., 
would violently expand into the outer flow, producing 
strong shock waves. Outward propagating shocks could 
heat the entire cooling flow provided the gas density has 
a steep radial dependence, n oc r~^ with s > 2 as in the 
models of Ciotti & Ostriker (2001). 

However, the (generally high redshift) EGRET sources 
exhibit large amplitude, random time variability. This 
variability is thought to occur because the radiation is 
beamed and Dopper boosted in a jet directed along our 
line of sight. This interpretation is supported by the as- 
sociation of gamma-dominant EGRET blazars with radio 
loud, flat spectrum central VLBI sources having large su- 
perluminal jet velocities. None of the EGRET sources are 
normal Seyfert galaxies, radio quiet QSOs or E galaxies 
except possibly Cen A (Hartman et al. 1999). Therefore, 
the assumptions that EGRET sources are isotropic gamma 
ray emitters and exist in all luminous elliptical galaxies are 
in doubt. A highly collimatcd beam of gamma radiation 
would heat the cooling flow much less efficiently. Nulsen 
& Fabian (2000) claim that the Compton temperature of 
more typical isotropic AGNs and QSOs is less than lO'' K, 
so cooling flow gas is more likely to be Compton cooled 
than heated. Nevertheless, it may not be possible to rule 
out isotropic bursts of gamma rays provided the duty cy- 
cle is sufficicintly short. If so, the Compton heating would 
be quite localized, decreasing as r~^ from the AGN. 

2.1.4. Shock Waves from Expanding Radio Plasma 

If the energy density in radio lobes is sufHciently large, 
they will expand supersonically, producing shocks that 
heat the cooling flow gas (Clarke, Harris, & Carilli 1997; 
Heinz, Reynolds, & Begelman 1998; Begelman & Cioffi 
1989; Rizza et al. 2000; Reynolds, Heinz, & Begelman 
2001; Soker, White, David & McNamara 2001). While 
such heating probably occurs at some level, recent Chan- 
dra and XMM observations show little evidence that the 
radio sources are actually heating the bulk of the gas or 
even gas just adjacent to the radio lobes (Fabian 2001; 
McNamara 2001). 

2.2. Recent Cooling Flow Models with Heating 



Loewenstein, Zweibel & Begelman (1991) discuss steady 
state cooling flows that are heated by cosmic rays in the 
center and by thermal conduction from larger radii. While 
their model is successful in producing central regions with 
dT/dr > 0, the central gas density gradient dn/dr is flat- 
ter than observed because the gas is both supported and 
heated by cosmic rays. They conclude that cosmic-ray 
heating cannot balance radiative cooling in intracluster 
cooling flows. 

Citing the difficulties encountered in producing thermal 
instabilities in nearly hydrostatic cooling flows (Malagoli, 
Rosner & Bodo 1987; Balbus & Soker 1989; Hattori & 
Habe 1990; MalagoH, Rosner & Fryxell 1990), Tabor & 
Binney (1993) assume that gas cannot cool in a spatially 
distributed fashion, and devise steady state flows in which 
a convective core heated from the center is attached to an 
outer cooling inflow. Unfortunately, in the convective core 
the temperature gradient is determined by the gravita- 
tional potential, {k/ finip) (dT/dr) = — (2/5)(d$/dr), and 
must always be negative, contrary to observation. Further- 
more, the net gas velocity within the convective core is ev- 
erywhere negative, causing all the gas shed from evolving 
stars to cool at the very center. After a Hubble time the 
mass of centrally cooled gas ~ 10^° Mq is unacceptably 
large. Binney & Tabor (1995) discuss time dependent, cen- 
trally heated models in which distributed mass dropout is 
not permitted, except perhaps near the core. These mod- 
els are heated by a large number of Type la supcrnovae at 
early times which can drive a galactic outflow throughout 
the evolution but which also enrich the hot gas with iron 
far in excess of observed abundances. Strong central su- 
pernova heating results in dT/dr < at every time in the 
flow evolution, unlike observed proflles. Their best-fltting 
models for NGC 4472, 4636 and 4649 agree with the ob- 
served gas density profiles only during a brief ~ 10^ year 
period just before a cooling catastrophe when the central 
gas density increases rapidly due to radiative losses. In 
additional models, Binney & Tabor include intermittent 
AGN heating within ~ 1 kpc triggered by gas inflow into 
the center. As SNIa heating subsides and the flow evolves 
toward a central cooling catastrophe, the AGN heating is 
activated but the central gas density still continues to in- 
crease, and, as Binney & Tabor suggest, this may result in 
distributed cooling and star formation in the central 1-2 
kpc. By the standards we set for success in our heated 
models, this flow would be regarded as a failure because 
cooling within a few kpcs would be visible as a multiphase 
component in XMM spectra. 

Tucker & David (1997) study time dependent models of 
cluster scale heated flows without cooling dropout or mass 
supply from stars. Gas near a central AGN is intermit- 
tently heated by relativistic electrons when gas flows into 
the origin. For these models, the temperature and entropy 
gradients are both negative while the density gradient is 
positive; such flows are both convectively and Rayleigh- 
Taylor unstable. When Tucker & David include thermal 
conduction, the outward heat flux results in stable solu- 
tions having cores of nearly flat density, but which are 
nearly isothermal, in disagreement with observations. In 
solutions with thermal conduction the mass flow M into 
the central regions is lower, and M{r) slowly varies for 
many Gyrs as a long-lived transient. 
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If gas in isolated elliptical galaxies is locally Comp- 
ton heated by powerful isotropic gamma rays {Lf,oi ^ 
1046 _ 1048 erg s-i; Tc ~ 10^ K) following mass accre- 
tion at the center, strong outward moving shocks can heat 
gas out to ^ 5 kpc and beyond (Ciotti & Ostriker 2001). 
The AGN-QSO continuum is only activated about ~ 10"'^ 
of the time, so only one in ~ 1000 of known elliptical galax- 
ies should be a gamma-dominant quasar at any time. For 
these models Ciotti & Ostriker assume a high Type la 
supernova rate and no supply of gas by secondary infall. 
Such cooling flows, which are less bound and have small 
inertial masses, can be significantly modified by the addi- 
tional AGN heating. In these models gas is driven out of 
the galactic core by shocks, so distributed cooling dropout 
is not necessary nor do large masses of gas cool at the 
galactic core; these are desirable attributes. As with previ- 
ous calculations with heating, however, dT/ dr is generally 
negative and the gas density is has a broad central core. 
The density and temperature profiles vary depending on 
the recent history of heating events during the previous 
sound-crossing time. In their preferred model 1 the cen- 
tral core in the X-ray surface brightness distribution 
varies in size from 1 to 40 kpc which is much larger than 
the ~ 0.6 kpc core in NGC 4649, an elliptical with an Lb 
similar to that used in the Compton heated models. Be- 
cause of radially propagating shock waves, the X-ray im- 
ages of these models are considerably more irregular than 
Tlx profiles observed in elliptical galaxies. As in the mod- 
els by Binney & Tabor (1995), the X -ray luminosity 
is expected to be large during those phases of the heat- 
ing cycle when the X-ray size is small and the X-ray 
image is more concentrated. But observations show just 
the opposite: and are positively correlated (Math- 
ews & Brighenti 1998). Finally, such large excursions in 
the Lx-rx plane following explosive heating events would 
be greatly reduced in more realistic cooling flows having 
larger masses of gas acquired from secondary infall, which 
extends far beyond the optical galaxy. 

3. CRITERIA FOR A SUCCESSFUL HEATED FLOW 

We seek computed solutions for heated cooling flows 
that satisfy three criteria. 

(1) After evolving for ~ 13 Gyrs, the hot gas density 
and temperature profiles should resemble those observed 
in massive E galaxies. For this purpose we adopt NGC 
4472 as a representative luminous E galaxy for which the 
density and temperature are well-determined from X-ray 
observations. 

(2) When our models are computed to t = Gyrs, as- 
sumed to be the current time, the rate of cooling (mass 
dropout) should not exceed cooling rates M observed in 
typical bright E galaxies. Unfortunately at the present 
time only a few galactic scale cooling flows have been ob- 
served with XMM and the indicated values of M based 
on spectral features differ quantitatively. Nevertheless, X- 
ray observers all agree that M is much lower than predic- 
tions of standard cooling flow models. From their XMM- 
EPIC observations of M87 Molendi & Pizzolato (2001) find 
M < 0.5 Mq yr-i while 40 ± 5 M© yr'^ and - 10 Mq 
yr~^ were predicted by Peres et al. (1998) and Fabian et 
al. (1984) respectively. More recently Matsushita et al. 
(2002) claim to have detected a total cooling rate in M87 



within 20 kpc of 0.8 Mq yr~^ while 8 Mq yr^^ was ex- 
pected in their traditional cooling fiow model inside the 
same radius, but the statistical likelihood of this solution 
is small. Recent XMM observations of the E galaxy NGC 
5044 fail to find spectroscopic evidence of cooling below 
~ 1 keV, implying an M consistent with zero (Buote et al. 
2002). For NGC 4636, which more closely resembles our 
fiducial galaxy NGC 4472, the XMM-RGS observations of 
Xu et al. (2001) can only provide upper limits to the total 
mass cooling rate. Xu et al. find a complete absence of 
emission lines from gas cooling below 0.5 keV. For an up- 
per limit they find M < 0.2 - 0.3 Mq yr"\ which is 3-5 
times less than cooling rates based on X-ray imaging, ^ 1 
Mq yr^^. However, FUSE observations of the OVI dou- 
blet by Bregman, Miller & Irwin (2001) indicate that gas 
in NGC 4636 is cooling through temperatures T ^ 3 x 10^ 
K at a rate M = 0.43 ± 0.06 Mq yr"^ In view of these 
inc;onsistencics in the ciirrent observations, we determine 
the success of a computed heated cooling flow in reduc- 
ing M by comparing with the cooling rate of an identical 
galaxy flow model without heating. Guided by the X-ray 
observations, the M of a heated flow should be reduced 
by about an order of magnitude, M ^ 0.1 Mq yr~^, to be 
compatible with XMM observations. 

(3) Finally, when our heated cooling flow models are com- 
puted to the current cosmic time, t ~ 13 Gyrs, the total 
mass of cooled gas within any radius must not exceed the 
combined mass of the stars and central black hole as de- 
termined from the stellar velocity dispersion. The black 
hole mass is - 2.6 x lO'' Mq for NGC 4472 (Magorrian et 
al. 1998). The stellar mass within about 200 pc is equal to 
the black hole mass. At a radius of 100 pc, within which 
much of the gas cools in many of our galactic flow models, 
the combined stellar and black hole mass in NGC 4472 (as 
described below) is M(r < 100 pc) ~ 2.9 x 10^ Mq. 

4. EQUATIONS 

Our reference galaxy NGC 4472 has a stellar mass 
M^t = 7.26 X 10^^ Mq, assuming M/Lb = 9.2 and a 
distance of 17 Mpc. If the X-ray luminosity of NGC 4472 
Lx{0-5 - 2 keV) — 4.54 x 10*^ erg s^^ is assumed to have 
been constant over tn = 13 Gyrs and there had been no 
heating, the total mass of cooled interstellar gas would be 
Mcooi ~ tnLx/{5kT/2iimp) » 4 x 10^° Mq. If there is no 
heating, the current rate of cooling in NGC 4472 is ^ 1.2 
Mq yr~^ (see below). 

For galactic scale cooling flows, the usual gas dynamical 
equations must include appropriate source and sink terms: 
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where e = 2>kT/2p,mp is the specific thermal energy, 
= 0.62 is the molecular weight and rup the proton mass. 
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The terms involving a = a* + a^n represent the addition 
of mass, momentum and energy to the hot gas by stellar 
mass loss and Type la supernova; usually <C a*. The 
stellar mass loss varies as a* = a*(t„)(t/f„)^^'^ s~"^ where 
a^{tn) = 4.7 X 10-20 s-i (e.g. Mathews 1989). 

The gravitational potential $ in the momentum equa- 
tion includes contributions from both stellar and dark 
mass. The de Vaucoulers stellar density p*(r) and mass 
M* (r) profiles are uniquely determined by the total mass 
M*t = 7.26 X 10" Mq and effective radius re = 8.57 kpc 
(e.g. Mellier & Mathez 1987). We assume a flat stellar 
core inside the break radius = 0.205 kpc (Faber et al. 
1997), which contains less than one percent of the total 
stellar mass. The NFW halo in NGC 4472 has a current 
virial mass Mh.vir = 4 x lO"'^"^ A4q based on its extended 
X-ray image (Brighenti & Mathews 1997). For the sim- 
ple f2 = 1 cosmology we consider with Hubble constant 
Ho = 50 km s~^ Mpc~^, the current dark halo density 
distribution is 

Ph ^ 200 1 

Pc 3 + c) - [c/ {I + c)] cx{l + cxY 

where pc = BH^/SttG is the critical density, c = 
4:98[Mh,vir/MQ]~°-'^'^^^ is the concentration, x = r/rh,vir, 
and rh,vir — (3Af;i,,,jr/47r200pc)^^^ (Navarro, Frenk & 
White 1996). While the large halo mass Mh,vir is appro- 
priate for the galaxy group from which NGC 4472 formed, 
our results are surprisingly insensitive to the (uncertain) 
value of Mh^vir- For example, when the same flow is com- 
puted with a halo mass reduced by ten, M^^yir = 4 x 10""^^ 
Mq, we typically find that the gas temperature and den- 
sity for r ~ 10—50 kpc are reduced by only about 15 and 45 
percent respectively. Since our assessment of the success of 
heated cooling flow models is weighted toward agreement 
with observations of NGC 4472 within the bright optical 
galaxy, r ^ 20 kpc, these relatively small changes in the 
flow models due to the uncertain value of Mh.vir have no 
effect on our conclusions. 

Guided by the relative sizes of the optical and X-ray 
images for E galaxies (Mathews & Brighenti 1998), we 
consider two types of evolutionary calculations for the hot 
galactic gas: "isolated" and cosmologically "accreting". 
Accreting models have large halos of X-ray gas acquired 
from cosmic flow into the galaxy group potential in which 
the central E galaxy formed. Alternatively, in our isolated 
E galaxy models all of the interstellar gas is supplied only 
by stellar mass loss and supernovae. In isolated galaxies 
the X-ray and optical images are of comparable size. In 
spite of our designation "isolated," (typically non-central) 
E galaxies with comparable X-ray and optical sizes prob- 
ably lack hot circumgalactic gas either because of tidal 
stripping or because they never received their full alloca- 
tion of cosmic gas. For simplicity, we assume that isolated 
E galaxies evolve from t = ti = 1 Gyr to the current 
time tn = 13 Gyrs with stellar and dark matter potentials 
fixed at their current values. In these models we begin 
the calculation by filling the galaxy with a very small gas 
density (n = 10~^ cm~^) at or near the virial temperature 
(T = 10^ K). Shortly after the calculation begins the mass 
that fiows from the stars overwhelms the initial gas and 
the final solution is independent of the initial gas. Accret- 
ing models are assumed to evolve from ti in a simple fiat 
universe, 0=1 with Hubble constant h = H/lOO = 0.5 



and Qh = 0.05, in which the dark halo grows from the 
inside out, beginning with a tophat cosmic perturbation 
(see Brighenti & Mathews 1999 for further details). In 
these models the continuing infall of cosmic gas creates 
hot gas distributions that extend far beyond the optical 
image. These fiows approximate E galaxies at the center 
of groups. 

The injection energy in Equation (3) is Sq = 3kTo/2p,mp 
where To — {a^T^ + asnTsn)/c(* is the mean temperature 
of gas expelled by stellar mass loss and supernovae. The 
supernova term 

asnTsn = 2.13xlO-«SNu(t) (j^^^) /j-^-'is"'' Ks'^ 

is proportional to the Type la rate in SNu units, SNu(t) = 
Q.{)Q{t/tny^ supernovae in 100 years per 10^"Lb. This 
value is consistent with the poorly known observed rate, 
~ 0.08±0.03 SNu from Cappellaro & Turato (2000). Since 
otsnTsn is generally much larger than a»T,, we simply as- 
sume T, = {jimp/k)al with a constant stellar velocity dis- 
persion (T, = 230 km s^^; the a,T* term has essentially no 
effect on our computed fiows. We evaluate agnTsn using 
Esn = 10^^ ergs; further details concerning the expression 
above for asnTsn can be found in Brighenti & Mathews 
(1999). Unless otherwise noted, a»(i), asn(i) and M** are 
the same for all niodc;ls and conduction and convection are 
zero, Fcond = ^conv = 0. If desired, we can include an ad- 
ditional heating source ff (r, t) in Equation (3) as described 
in more detail below. 

The dimensionless coefficient q in the sink term in Equa- 
tion (1) can be used to remove hot gas from the fiow 
over some large volume or to remove gas that has already 
cooled. In either case —qp/tcooi oc — gp^, where tcooi = 
5mpkT/2ppA is the time for gas to cool locally by radia- 
tive losses at constant pressure. For simplicity we evaluate 
the cooling coefficient A(T) using the results of Sutherland 
& Dopita (1993) assuming solar abundance. The cooling 
rate and the gas density within ^ r,, are slightly underes- 
timated in central galaxies where zp^ ^ zp^.Q- 

Ideally, if the gas in heated cooling fiows does not cool, 
there would be no reason for this ad hoc sink term. How- 
ever, in the heated cooling flows models discussed below we 
find that gas cools to low temperatures in localized regions 
(drips) where radial compression (V • u < 0) is sustained. 
To remove the accumulated cooled gas from the flow in a 
numerically satisfactory manner in spherical ID flows we 
let 

q{T) = qcooi eM-T/Tcooif (4) 

where qcooi = 2 and Tcooi = 5 x 10^ K; we shall discuss the 
sensitivity of the flow solutions to these two adjustable pa- 
rameters. When T » Tcooi, as usually occurs, the cooling 
dropout term is negligible and has no influence on the flow. 
We use Equation (4) mainly to locate where gas is cool- 
ing; physically, it could represent the formation of stars 
or dense clouds that orbit near their radius of formation. 
Flows calculated with q = Q and a minimum tempera- 
ture Tmin = 10"* K are essentially identical to those using 
Equation (4); in this case gas that cools to Tmm = 10^ K 
freefalls into the center of the galaxy. While our results are 
insensitive to the precise mathematical form of Equation 
(4), in our 2D calculations we prefer to remove cooled gas 
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with a step function: 

where Tcooi = 5 x 10'* K. This step function with its lower 
Tcooi and sudden onset better illustrates the geometrical 
complexity of of the cold gas distribution in 2D flows. 

For ID and 2D models we use a modified version of the 
ZEUS-2D code (Stone & Norman 1992). In ID there are 
300 zones within Vmax = 500 kpc that slowly increase in 
size from the smallest zone at 50 pc. To capture the com- 
plex small scale structures of out 2D calculations, we use 
a cylindrical grid with 600 zones along each axis oTit to 
Zmax = Rmax = 190 kpc; the central zone has dimensions 
^1 = i?i = 100 pc. 2D solutions with 120 x 120 or 300 x 300 
zones (and z\= R\ = 100 pc) are qualitatively similar. 

The ZEUS code has been widely tested and is known to 
accurately compute a variety of standard hydrodynamical 
problems (Stone & Norman 1992). Our detailed source 
term modifications of ZEUS (Equations 1-3) have been 
checked by comparing our galactic flow models with iden- 
tical calculations using the independent cooling flow hy- 
drocode developed by Annibale D'Ercolc in Bologna (e.g. 
D'Ercole & Ciotti et al. 1998). Finally, we have suc- 
cessfully checked the 2D ZEUS code against analytic ex- 
pressions for nonsteady inflow in hypothetical galaxies of 
uniform density (Appendix of Brighenti & Mathews 1996). 

In the following we discuss a large number of ID and 
2D flows, both heated and unheated. Readers who wish 
to review a representative sample of each type of model 
may elect to skip §6.3 and 6.4. However, the ID and 2D 
flows described in §6.3 and 6.4 are essential in establish- 
ing the generality of our conclusion that heating cannot 
substantially reduce the global cooling rate in the flows. 

5. COOLING FLOWS WITHOUT HEATING 

For comparison with heated models it is useful to first 
calculate traditional cooling flows without heating or other 
complications, i.e. H{r, t) = Fcond = Fconv = 0. 

5.1. No Cooling Dropout, q = 

The left column of Figure 1 shows the gas density, tem- 
perature and velocity at t„ = 13 Gyrs for single phase 
g = flows in both isolated and accreting galaxies. It is 
seen that the NGC 4472 data can be matched reasonably 
well with the isolated solution for 1 < r < 10 kpc, but 
the temperature is low and has the wrong gradient, i.e. 
dT/dr < 0. In r > 10 kpc the gas density is also too 
low in the isolated galaxy. The outer density profiles are 
better fit with the accreting solution (dashed lines), indi- 
cating that not all of the hot gas in massive E galaxies 
comes from stellar mass loss (Brighenti & Mathews 1999). 
Of particular interest for the q = accreting solution in 
Figure 1 is its agreement with the positive temperature 
gradient observed beyond about 10 kpc. In all well-studied 
galactic cooling flows the gas temperature rises from a 
minimum at the origin to a maximum near ~ (3 — 4)re 
where the temperature is about 30 percent larger, then 
slowly decreases at larger radii beyond the optical image of 
the central galaxy (Brighenti & Mathews 1997). High gas 
temperatures at large radii are produced in the accretion 
shock when the gas first entered the massive dark halo of 



the galaxy group surrounding the luminous elliptical. The 
post-sliock temperature of the accreting gas is close to the 
virial temperature of the entire galaxy group. The pos- 
itive temperature gradient within several Ve is produced 
as the temperature approaches the somewhat smaller ^ 1 
keV virial temperature of the central galaxy. This temper- 
ature profile is consistent with that of NGC 4472 shown 
in Figure 1. Finally, we note that the temperature T{r) 
plotted in Figure 1 and in subsequent figures is shown as 
a function of physical radius r, not the slightly different 
emission-averaged temperature along the line of sight at 
the projected radius. Although this latter temperature 
profile would be more appropriate to compare with the 
NGC 4472 observations, T as a function of r is necessary 
to facilitate comparisons with ID heated cooling flows (§6 
below). 

The computed density n(r) in 0.1 ^ r ^ 1 kpc for both 

q = solutions exceeds the observations; indeed, the hy- 
pothesis of distributed cooling (mass dropout) g > was 
first devised to correct this discrepancy (e.g. Thomas 1986; 
Vedder, Trester & Canizares 1988; Sarazin & Ashe 1989). 
However, the most obvious and fatal problem with 5 = 
models is that 4 x 10*° Mq of cooled gas has accumu- 
lated in the central computational zone by time t„. As a 
consequence of this mass, the hot gas density inside 100 
pc rises sharply as seen in Figure 1. This central mass is 
~ 20 times larger than the black hole mass observed in 
NGC 4472 and causes the gas temperature to peak unre- 
alistically in r ^ 5 kpc, unlike any known bright elliptical 
galaxy. 

As a historical point, in earlier steady state galactic cool- 
ing fiow models the global temperature profiles were sensi- 
tive to the outer boundary conditions. White and Cheva- 
lier (1984) required that the gas density vanish with the 
flow velocity at some large stagnation radius rgtag- This 
constraint causes the temperature to rise with radius and 
become infinite at rgtag^ assuming the gas pressure there 
is finite. This sensitivity of steady state flows to the outer 
boundary conditions was recognized by Vedder, Trester & 
Canizares (1988) whose solutions - unconstrained by the 
density at rgtag reveal a wide variety of global temper- 
ature profiles. They realized that steady state boundary 
conditions imposed at Vatag are unrealistic since the flow 
time to the origin typically exceeds the Hubble time. 

In summary, isolated flows without cooling dropout have 
too little density beyond ^ 10 kpc and have negative 
rather than positive gas temperature gradients. Accret- 
ing flows are in good agreement with the extended hot gas 
density structure and the positive temperature gradient 
within several r^- But both g = solutions fail because 
the central mass of cooled gas is too large. 

5.2. With Uniform Cooling Dropout, q = 1 

Also shown in Figure 1 are traditional cooling flows with 
distributed (multiphase) cooling dropout at every radius, 
g = 1 (e.g. Sarazin & Ashe 1989). Both the density and 
temperature are reasonably well matched, except in the 
central 1-2 kpc. The temperature derivative is positive 
for the accreting model in 5 < r < 50 kpc in agreement 
with observations. The temperature is slightly lower than 
observed in NGC 4472, but this could be corrected by 
small adjustments in the parameters of the calculation 
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(e.g. Brighenti & Mathews 1999). The NGC 4472 ob- 
servations are included in the Figures only for reference 
and we have not attempted to fit them exactly. 

The central density in both q = 1 solutions is still 
too peaked, but a flat density core could be produced by 
adding magnetic or cosmic ray pressures not included in 
Figure 1. The mass that cools at the origin is only 2.1 x 10^ 
Mq and 2.9 x 10^ Mq for the isolated and accreting q = 1 
flows respectively, but the total mass of cooled gas is very 
similar to the q ^ flows. Additional details on models of 
this sort can be found in Brighenti & Mathews (2000). 

While the q= 1 dropout solutions are generally success- 
ful, the distributed dropout hypothesis has been criticized 
for several reasons: (1) there is as yet no XMM-Chandra 
spectral evidence for cooling to T = 0, (2) linear perturba- 
tions do not result in localized cooling, (3) large masses of 
colder, cooled gas have not been observed, and (4) contin- 
uous star formation requires deviations from the Salpeter 
IMF. 

6. HEATED COOLING FLOWS 

6.1. When Heating Locally Balances Cooling 

Some authors have noted that the observed mass of hot 
gas in largo E galaxies is roughly equal to the total mass 
ejected from stars, Mgas ~ M*(a*(t„)t„ ~ 10^" Mq, and 
have concluded from this that there is no need for mass 
dropout as long as some hypothetical heating source bal- 
ances radiative losses at every radius. In the left column 
of panels in Figure 2 we show idealized solutions of this 
sort without distributed cooling (g = 0) and assuming 
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so that radiative cooling is effectively zero. In these solu- 
tions u « and no gas cools at any radius. However, the 
gas density is very large at time tn = 13 Gyrs and much 
too centrally peaked [p oc p*) due to the absence of an 
inward flow. As the mass of interstellar gas accumulates, 
a larger fraction of the pressure gradient required for virial 
equilibrium is taken by the density gradient so the tem- 
perature and its gradient are somewhat low. We have also 
computed flows for which the dissipation of stellar mass 
loss ap*eo,* is the only heating, i.e. 
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pH{r, t) = — 5" - ap, 
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where eo,sn = {?>k /2iimp)(asnTsn) I oi*, with results very 
similar to those shown in Figure 2. While this type of glob- 
ally heated flow might seem to be the most perfect possible 
solution to the mass dropout problem, it fails miserably to 
fit the observed data. 

6.2. A Dripping Flow with Distributed Heating 

Before we discuss results of a large number of flow cal- 
culations with various assumed heating terms H{r,t), it 
is usehil to consider in detail a specific but representative 
fiow with nonzero heating. For this purpose we adopt a 
Gaussian profile: 

H{r,t) = GheM-ir/rhf] (8) 
with Gh = I erg gm~^ s~^ and Vh = 5 kpc. The heating 
is instantaneously applied whenever the gas velocity into 



the central zone is negative, otherwise H{r,t) = 0. The 
heating feedback tends to inhibit gas from flowing into 
the center. This heating scenario is clearly a simpliflca- 
tion, but if additional features were included - time delay 
for feedback, slower rise or decay time for H(r,t), slower 
propagation of heating through the gas, etc. - these details 
would obscure but not significantly alter our main results. 
As noted earlier, we do not attempt here to justify H{r, t) 
in terms of any particular physical model for heating, our 
approach is purely phenomenological. 

In the right column of Figure 2 we show the flow 
variables for both isolated and accreting models at time 
tn = 13 Gyrs as influenced by the Gaussian heating of 
Equation (8). At earlier times, t < tn, the flow velocity 
into the central zone is generally negative, so the heating 
H{r,t) is almost continuous. The gas density is globally 
lower than observed in NGC 4472, but the shape is rea- 
sonably correct, particularly for the accreting fiow. How- 
ever, the gas temperature in the accretion flow is too large 
for r > 5 kpc due to the central heating and the slope 
dT/dr < 0, is wrong. 

The striking new feature in this Gaussian-heated flow 
is a nearly stationary "drip" wave at r « 7 kpc visible in 
Figure 2 where the gas density has a sharp maximum and 
the temperature drops suddenly to very low values. (This 
is why we do not attempt to plot the average temperature 
along the line of sight.) This type of concentrated drip 
cooling is a standard attribute of ID spherically symmet- 
ric heated models. The drip in the right column of panels 
in Figure 2 occurred because the heated gas within 7 kpc 
is slowly flowing out, while the negative velocity in r > 7 
kpc behaves much like a normal cooling flow. Just at r ~ 7 
kpc the gas is locally compressed, V ■ u < 0, generating 
runaway radiative cooling. 

The global velocity field has produced a spontaneous 
nonlinear thermal instability far from the center of the 
fiow. We note that "drips" were first discovered in un- 
heated cooling flows (Mathews 1997) where the differential 
(negative) flow velocity can occasionally produce a local 
compression. If cooled gas is not removed from the drip, 
the mass of locally cooled gas increases until, in ID cool- 
ing flows, it goes into freefall and collapses supersonically 
toward the center of the flow, producing reflected shocks 
from the origin and altering the global density and tem- 
perature profiles in very unrealistic ways. Our drips move 
slower or even become approximately stationary since we 
remove cooled gas inside the drip by adopting q > (Equa- 
tion 4). 

Drips may have been missed or misinterpreted in previ- 
ous discussions of heated cooling fiows. They are particu- 
larly easy to capture in Lagrangian calculations (Mathews 
1997), but can be missed in Eulerian calculations particu- 
larly if the grid is too coarse. One-dimensional drip waves, 
such as those shown in Figure 2, are clearly unphysical 
since they are unstable both to Rayleigh- Taylor disrup- 
tion (Vp opposes g = — V$) and to convective motions 
(negative entropy gradient inside the drip structure). The 
true evolution of drips can be understood only with 2D or 
3D calculations, discussed below, but the incidence of spa- 
tially distributed cooling, represented in ID with a drip, as 
well as the total mass and approximate location of cooled 
gas, are determined reasonably well from ID calculations. 
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We find for the isolated ID galaxy (riglit column of Fig- 
ure 2) that 4.0 X 10^ Mq has cooled for r < 100 pc and an 
additional 6.1 x 10*^ Mr.^ has cooled in 0.1 < r < 9.5 kpc 
due to inward propagating drips during t < 3 Gyr. Figure 
3 shows the mass distribution of cooled gas Mcoo;(^) fo'^ 
the isolated heated flow at two times, 3 and 13 Gyrs, as- 
suming that the cooled gas remains at its radius of origin. 
We cannot fully trust the astronomical accuracy of our 
flow solutions during the first 1-2 Gyrs of the calculation 
{t < 3 Gyrs). Since the stellar mass loss is particularly 
copious at these early times, it is possible that relatively 
large masses of cooled gas combined and formed stars with 
a normal IMF. By the current time t„ these stars would 
be almost as old at the original population and would con- 
tribute to the luminous stellar mass within several kpc. So 
the large mass 4.0 x 10^ M© that flowed into the central 
grid zone before t = 3 Gyrs in this calculation, which is 
larger than any known galactic black hole, cannot with 
confidence be used to disqualify this model. For this rea- 
son we distinguish in Figure 3 the mass that cooled before 
and after t = 3 Gyrs. 

After time t = 3 Gyrs almost all of the cooling gas 
in the accretion or isolated flows is formed in the nearly 
stationary drip at 7 kpc. The concentration of cooling 
gas at this radius is an artifact of spherical symmetry; in 
an otherwise identical 2D flow (discussed below) we find 
more distributed cooling in r < 10 kpc, but the total mass 
of cooled gas is unchanged. The cooling rate during the 
last 3 Gyrs is ~ 1 Mq yr~^. This nearly continuous drip 
cooling would produce a multiphase spectrum in apparent 
disagreement with X-ray spectra observed with XMM in 
M87 and other similar large E galaxies. Ironically, there- 
fore, the heating process that has been invoked to eliminate 
distributed cooling in cooling flows has instead stimulated 
widespread cooling just beyond the region of intense central 
heating. 

Another undesirable characteristic of heated cooling 
flows is the negative temperature gradient within several 
tens of kpc from the galactic center, as indicated in the 
right column of Figure 2 for r > 5 kpc. This variation is 
opposite to the thermal profile observed in all luminous, 
group-dominant elliptical galaxies. 

Thermal conduction is likely to be important within the 
sharp ID drip structure where the temperature gradients 
become large. Convection must also occur since the en- 
tropy gradient is locally negative. Both processes are ex- 
pected to broaden the drip structure. But the question of 
most interest is whether the cooling dropout can be slowed 
or completely stopped by including conduction or convec- 
tion. We show in the Appendix that the mass dropout is 
not significantly reduced by either conduction or convec- 
tion. 

6.2.1. The Same Flow in Two Dimensions 

Spherically symmetric ID drip waves are inherently un- 
realistic because of local density and entropy inversions. 
We now describe the more realistic two dimensional evo- 
lution of the same Gaussian heated flow for the isolated 
galaxy assuming Fcond — Fconv — 0. We use cylindrical 
coordinates since this geometry does not favor the spheri- 
cal character of the initial drip evolution, i.e. if distributed 
drip cooling occurs in 2D cylindrical flow, it would cer- 



tainly be expected in 2D spherical coordinates. Figure 
4 shows the density and velocity structure at tn for the 
2D version of the isolated galaxy flow in the right col- 
umn of Figure 2 heated as described by Equation (8) with 
Gfi = 1 erg gm~^ and rh = 5 kpc. Since the flow is 
heated from below, subsonic turbulence is generated near 
the center, r ^ 25 kpc. 

The most interesting result of the 2D fiow is that strong 
spatially distributed mass dropout still occurs. Also shown 
in Figure 4 are contours of the cooled gas density gen- 
erated with the step function q{T) in Equation (5) with 
Tcool = 5 X 10^ K, again assuming that the cooled gas 
remains where it first appears in the flow. Because of 
Rayleigh- Taylor and convective instabilities, the distribu- 
tion of cooled gas is no longer confined to a narrow drip 
region as in the corresponding ID flow (Figure 2), but cool- 
ing still occurs in about the same region. The tendency 
for additional cooling to occur along the R and z axes is a 
numerical effect: when convective eddies flow toward the 
R = OT z — boundaries, the reflection boundary condi- 
tions result in local compressions and subsequent cooling. 
However, widespread cooling clearly occurs throughout the 
2D flow within radius r < 10 kpc when convective eddies 
encounter each other, producing local nonlinear compres- 
sions. (A similar localized compressive cooling far from 
the galactic center also occurred in the calculations of gas 
flows in SO galaxies by D'Ercole & Ciotti 1998). By time 
tn the total amount of cooled gas Mcooi = 3.8 x 10^° Mq is 
comparable with that found in the corresponding ID flow. 
Azimuthally averaged hot and cooled gas densities, n(r) 
and ncooi{r), for this 2D flow are also shown in Figure 4. 

The maximum radial extent of the cooled gas in 2D flows 
increases as the critical dropout temperature in Equation 
(5) Tcool is increased above 5 x 10"* K. This can be under- 
stood in terms of our ID solutions in which drips typically 
begin at larger radii and slowly migrate inward as their 
amplitude and mass increase. The temperature minimum 
in the drip structure deepens as the drip moves inward. 
Consequently, increasing Tcool in 2D flows systematically 
removes gas at somewhat larger radii. Because of this 
sensitivity to the parameter Tcool and to the grid resolu- 
tion, we cannot be confident of the precise distribution of 
cooled gas shown in Figure 4, but the reality of spatially 
distributed cooling (dropout) in these heated flows is ro- 
bust. The Gaussian heated cooling flow shown in Figure 
4 would therefore have a multiphase spectrum (and tem- 
perature proflle) unlike those observed with XMM. 

In view of the overall similarity of ID and 2D calcula- 
tions of heated cooling flows such as those in Figures 2 
and 4, we assume that if cooling drips occur in ID, dis- 
tributed cooling will also appear in a more realistic 2D 
flows having the same parameters. We have not proven 
this assumption, but it seems to hold true for the cases we 
have checked. As a result, we can explore a wide variety of 
heated cooling flows in ID and assume that a comparable 
mass dropout would also appear in similar 2D calculations. 
The ID heated flows are not only quicker to compute, they 
show more clearly where and why the initial cooling oc- 
curs; the corresponding 2D flows reveal more realistically 
the geometrical complexities of the subsequent cooling. 

To illustrate the similarity of global cooling rates in sim- 
ilar ID and 2D fiows, we compare in Figure 5 the total 
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mass of cooled s function of time Mcooi{t) for the 

ID and 2D Gaussian heated models. The agreement is 
very good. It is also remarkable that the total mass of 
gas that cools and the mean cooling rate during the last 
few Gyrs are insensitive to the presence of conduction or 
convection (see Appendix), the cooling dropout threshold 
Tcooi, and the computational geometry, ID or 2D. In 2D 
the distributed cooling is intermittent, occurring in bursts 
alternating with quiescent periods of several Gyrs without 
dropout. 

Cooling intermittency in 2D calculations can be ex- 
plained by a quasi- cyclical variation of the density profile. 
After each cooling episode, the hot gas density is lowered 
and turbulence-driven compressions are insufficient to pro- 
duce local thermal instabilities. Later, after 1 — 2 Gyrs 
when the hot gas density has increased due to stellar mass 
loss and inflow, post-compression radiative losses become 
more efficient in triggering cooling dropout. In Figure 6 we 
show the evolution of the hot gas density profile during a 
cooling episode for the 2D ffow shown at time t„ in Figure 
4. During the cycle the ROSAT band X-ray luminosity 
varies from = 1.42 x 10^^ ergs s~^ at time t = 8 Gyrs 
to a maximum of = 2.66 x 10'*^ ergs at t — 9.8 
Gyrs and then back to = 8.76 x lO^*' ergs s~^ at t = 11 
Gyrs. Unfortunately, this variation of by about 30 is 
unlikely to be the source of the large scatter observed in 
the (La;,L_B)-plane since cooling is not supported by XMM 
observations and the large excursions in the density (and 
X-ray surface brightness) profile shown in Figure 6 are not 
typically observed (e.g. Trinchieri, Fabbiano & Canizares 
1986). 

6.3. More Flows with Distributed Heating 

In addition to the representative heated cooling flow 
that we just described, we have computed a large number 
of additional flows for both isolated and accreting galaxies 
with various heating scenarios H{r,t). Most of these are 
ID flows, but we assume that the corresponding 2D solu- 
tions are qualitatively similar. Our conclusion is that none 
of these flows is acceptable. Heated cooling flows typically 
fail one or (usually) more critical tests: (1) significant on- 
going cooling occurs within ~ 10 — 20 kpc which would be 
observed as a multiphase component in XMM spectra, (2) 
the temperature gradient in 1-10 kpc is negative or peaked 
in heated accretion flows, not positive as observed, (3) the 
density at 5-20 kpc is often too low in heated flows due to 
the redistribution of heated gas. While it is necessary and 
important to carefully document the shortcomings of each 
type of heated galactic flow, readers with a low tolerance 
for computational detail and who accept these general con- 
clusions may wish to skip to §7. 

We consider a wide variety of heating proflles beginning 
with the diffuse Gaussian profile in Equation (8) but with 
different parameters Gh and r^. A more concentrated in- 
verse square heating is described with 

H{r,t)=ChTcr-^ (9) 

where the time dependence arises because heating is 
turned off if there is no mass flow into the central grid zone. 
This one-parameter {ChTc) heating function is designed 
to simulate inverse Compton heating if Tc ^ 10^ 
K, but it can represent any form of unspecified concen- 
trated heating in the galactic core. To simulate Compton 



heating when Tq is closer to T^n we have also performed 
calculations with 

H{r,t)=Chr-\Tc-T). (10) 
In this case Ch = {4:aT/-i'meC^){3k/2ij,mp){Lagn/4:''^) erg 
cm^ g~i K"-'^ is an approximation to isotropic Comp- 
ton heating using the Thompson crossection. 

We also consider a heating term proportional to the stel- 
lar density. This type of heating was performed by sim- 
ply increasing the current Type la supernova rate SNu(t„) 
above the value 0.06 SNii that agrees with currently ob- 
served iron abundances (see Figure Al in Boute 2000). 
This type of heating Hp cx p*{r) is continuous, indepen- 
dent of inflow into the central grid, but H can also share 
the Type la time dependence, SNu oc t~^, so that SNu(t„) 
and s are the two adjustable parameters for this heating 
scenario. While solutions with s = 1 and SNu(t„) > 0.06 
clearly produce too much iron, our intention is to ex- 
plore the possibility of some type of hypothetical star- 
related or extended heating possibly unrelated to super- 
novae. Finally, we computed flows in which mass dropout 
is accompanied with simultaneous local heating feedback, 
pH{r,t) oc p = q{r)p/tcooi- 

6.3.1. Additional Gaussian Heated Profiles 

In the flrst two columns of Figure 7 we illustrate 
Gaussian-heated flows with r/j = 5 kpc as before but 
with larger and smaller heating coefficients Gh- Solid and 
dashed lines represent isolated and accretion flows respec- 
tively. For Gfi = 0.5, half as large as that in Figure 2 
(right column), the heating is insufficient to induce sig- 
nificant large scale drip cooling so most of the gas cools 
within the central 100 pc. For the isolated flow this mass 
is M,,ooi{r < 100 pc), is 1.0 x 10^° Mq and 3.0 x 10^° Mq 
at 3 and 13 Gyrs and for the accretion model the masses 
cooled within 100 pc are 0.74 x 10^° Mq at 3 Gyrs and 
3.9 X 10^° Mq at 13 Gyrs. These high central masses 
are clearly unsatisfactory. For the isolated flow the hot 
gas density is generally too low and for the accretion flow 
the gas density gradient is too steep, although the density 
profiles are somewhat time-variable. For both flows the 
gas temperature becomes large near the center due to the 
cooled mass there; as a result dT/dr is negative within 5 
or 10 kpc, not positive as observed. 

For Gh = 'i', twice as large as that in Figure 2 (right 
column), the larger heating causes the gas density to be- 
come much lower than the NGC 4472 density, but a large 
amount of gas still cools within 100 pc in both isolated 
and accretion flows, about 2 x lO'' M©, comparable to the 
black hole mass in NGC 4472. The gas temperature in the 
inner galaxy r < 10 kpc is far larger than observed and 
dT/dr < for r > 5 kpc. Continuous drip cooling near 
10-12 kpc would produce observable multiphase spectra 
for either accretion or isolated flows. In a more realistic 
2D version of these = 2 calculations we expect that 
this cooling would be distributed within ~ 15 kpc, not 
concentrated in the ID drip as shown in Figure 7 (second 
column). Clearly, the Gh parameter must be fine-tuned 
even to achieve a poor fit to the observations. 

Next we vary the spatial scale rh of the heating, then 
increase Gh until the heating begins to affect the flow. The 
accretion and isolated flows have similar problems as be- 
fore. For Gh = 0.5 and rh = 20 kpc (third column of Fig- 
ure 7), subsonic outflow occurs at r > 1 kpc, but the gas 
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interior to this radius flows to the origin most of the time, 
and the central mass is still much too large, 6 — 10 x 10^ 
Mq. Drip cooling occurs at r ~ 60 — 100 kpc, which would 
produce a multiphase spectrum. The gas temperature is 
too high beyond 5 kpc for both the isolated and the ac- 
cretion flow. For both types of flow the density is too low 
at all radii r > 1 kpc. Increasing Gh further lowers the 
density even more without appreciably reducing the cen- 
tral cooled mass. As expected, values of Gh less than 0.5 
(keeping = 20 kpc) result in strong central cooling and 
even larger central masses. 

When the Gaussian heating is more concentrated, r/j = 
1 kpc then Gh must be 20 to influence the flow (right col- 
umn of Figure 7). For these heating parameters, a central 
subsonic outflow is present at t„ within ~ r^, resulting 
in sustained drip cooling at r ~ 1.5 — 2 kpc. Interior to 
the drip the gas temperature is much too high, T ~ 2 — 4 
keV. This same pattern occurs in both accretion and iso- 
lated flows. While 3.7 x 10^° Mq cools in 1.5 < r < 2 
kpc, very little gas cools within 100 pc in these models. 
The central density is flattened with a rather good flt to 
NGC 4472 at all relevant radii r ^ 10 kpc, but is about 
three times too large for 2 < r < 6 kpc in both accretion 
and isolated flows. Overall, models with more centrally 
concentrated heating are in better agreement with obser- 
vations on large scales, but the deposition of extremely 
large masses of cooled gas within a few kpc and the as- 
sociated multiphase XMM spectra are the same undesired 
consequences encountered in other flows with hypothetical 
heating. 

6.3.2. Concentrated Heating with H rx 

In this series of flows the heating is described by Equa- 
tion (9). For cither isolated or accretion flows with 
ChTc — 1.6 X lO^'' erg cm^ g^^ s^^ (where for example 
Lagn « 10*^ erg s-i and Tc = 10^ K) the density gradi- 
ent at time tn is too steep (as shown in the left column of 
Figure 8). About 4 x 10^° Mq cools in r < 0.5 kpc, which 
results in a high temperature core within 1 kpc that is 
rarely if ever observed. If ChTc is increased to 1.6 x 10^'' 
to reduce the mass that cools near the origin, the central 
gas density is too low and the flat density core is unaccept- 
ably large (right column of Figure 8). The central temper- 
ature is also too high and multiphase cooling is expected 
at about 1.5-3 kpc. As before, results are similar for ei- 
ther accretion or isolated models. If thermal conduction is 
added to distribute the thermal energy to larger radii, nei- 
ther the density or temperature profiles are substantially 
improved and drip cooling still occurs but at a somewhat 
larger radius. In summary, as the cooling parameter ChTc 
increases by a factor of 10, the solutions shift from ones in 
which too much gas accumulates near the origin to flows in 
which the gas density is unrealistically lowered by strong 
heating. It is possible that an intermediate value of ChTc 
would produce a satisfactory density proflle, but the tem- 
perature would still be too high and multiphase gas near 
the center seems unavoidable. The flne-tuning required, 
moreover, seems very implausible. 

6.3.3. Distributed Heating with Hp oc p* 

It has often been suggested that some source of stel- 
lar heating, such as with Type la supernovae, could bal- 



ance radiative losses in cooling flows, arresting the cool- 
ing process. Therefore, we briefly consider ID isolated 
and accretion flows in which H{r)p oc p*{r). This can 
be accomplished by artificially increasing the supernova 
rate SNu(i) = SNu(t„)(t/t„)-* by increasing SNu(t„) or 
adjusting its time dependence. We describe several solu- 
tions for isolated and accreting flows in which SNu(t„) is 
increased while keeping s = 1. The transition from cool- 
ing concentrated near the galactic center to galactic winds 
is very sensitive to SNu(t„). Three isolated solutions are 
shown in the left column of Figure 9. For SNu(f„) = 0.20 
there is a cooling inflow and Mcool = 2.3 x lO^'' Mq cools 
within the central 100 pc; with SNu(t„) = 0.23 there is 
a partial outflow with cooling dropout near 3 kpc; with 
SNu(t„,) = 0.25 there is a complete wind and no gas cools 
at any radius. The accreting flows shown in the right 
hand column of Figure 9 exhibit the same sensitivity to 
SNu(i„) with inflow at SNu(i„) = 0.25, partial outflow 
at SNu(t„) = 0.28 and a fuU wind at SNu(t„) = 0.30. 
For either type of flow, the temperature gradient is neg- 
ative everywhere and the hot gas density suddenly drops 
to unacceptably low values as SNu(t„) exceeds a sharp 
threshold for wind solutions. In addition to its inability 
to reproduce the observed temperature and density pro- 
files, this type of heating with H{r) cx cannot be due 
to Type la supernova - if it were, the iron abundance 
Zpe would be extremely high; for example in the accreting 
SNu(f„) = 0.25 solution zpe would be more than 6 times 
solar within about 60 kpc. 

6.4. Additional 2D Heated Cooling Flows 

While the initiation of nonlinear thermal instabilities 
that appear in heated cooling fiows are clearly revealed 
in the ID flows, the later evolution of these instabilities 
is inherently multidimensional because of Rayleigh- Taylor 
and convective instabilities. For this reason we discuss 2D 
versions of several heated cooling flows already described 
with spherically symmetric models. For simplicity, we con- 
sider only 2D isolated flows. To implement the feedback in 
2D, heating is applied only when the mass flux at radius 
r = 300 pc is negative, otherwise we set H = (models 
with continuous heating give qualitatively similar results). 
Cooled gas is removed using the step function procedure 
(Equation 5) with T^ooi = 5 x 10^ K. 

The left column of panels in Figure 10 shows the 2D 
version of the Gaussian-heated flow with Gh = 20 and 
r/i = 1 kpc corresponding to the solid line ID solution in 
the right column of Figure 7. The upper panel shows the 
hot gas density and velocity and the second panel shows 
contours of the density of cooled gas, both at time i„. The 
lower two panels show the azimuthally averaged hot and 
cooled gas densities respectively. While most of the cool- 
ing along the vertical {z) axis is a numerical artifact, as 
explained earlier, there is genuine small-scale cooling away 
from both axes at r ^ 1 — 2 kpc. This cooling is visible 
as a peak in the azimuthally averaged cooled gas density 
at r 1.5 kpc in the bottom panel. Cooling in this re- 
gion is expected from the ID isolated flow illustrated in 
the right column of Figure 7. The total mass of cooled 
gas is 3.0 X 10"'^'^ Mq. The azimuthally averaged hot gas 
density distribution at time i„ shown in Figure 10 is much 
too steep relative to NGC 4472 and is steeper than the 
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corresponding density profile in Figure 7. We believe that 
this latter difference can be explained by the quasi-cyclic 
variation in the density profile associated with episodic 
cooling as shown in Figure 6; i.e. the ID and 2D flows 
are captured at different phases of this cycle at time t„. 
Like the corresponding ID flow, the 2D gas temperature 
(not shown) is also somewhat too high in the central kpc, 
T ~ 1.5 keV. This 2D heated cooling flow suffers from 
the usual astronomical problems: the mass of cooled gas 
within a few kpc is very large, radiatively cooling gas in 
this region should produce a multiphase XMM spectrum, 
and the variable excursions in the hot gas density profile 
are larger than the narrow range of density profiles consis- 
tent with X-ray observations (e.g. Fabbiano, Trichieri & 
Canizarcs 1986). 

The central and right hand columns of Figure 10 show 
2D flows at time f„ corresponding to the isolated ID flows 
illustrated in Figure 8 for which the feedback heating is 
Compton-like, H cx r^^. The 2D solutions confirm our 
previous conclusions for these flows: for weak heating, 
ChTc = 1.6 X 10^3^ aU the gas (~ 4 x 10^° Mq) cools 
at the center of the flow (central column Figure 10) and 
the gas density profile is too steep. If the heating is 
stronger, ChTc — 1.6 x 10'''', (right column Figure 10) 
widely distributed cooling occurs that is not seen in XMM 
spectra. Even with this larger ChTc, an excessively large 
mass 2.2 x 10^° Mq still cools throughout the galaxy. The 
central gas density is unrealistically low within about 2 
kpc and too high beyond. The temperature is too high 
(T ~ 1.5 - 2 X 10'^ K) inside ~ 10 kpc, similar to the 
ID solutions in Figure 8. Extended periods during which 
there is very little radiative cooling at large r are punctu- 
ated by short episodes of intense cooling, but the average 
cooling rate history is similar to the other models shown in 
Figure 5. Nevertheless, it is obvious that neither of these 
models with H cx is acceptable, nor is it plausible that 
any value of ChTc could give acceptable results. 

Finally, we have computed ID and 2D versions of cool- 
ing flows in which the cooled gas is rapidly heated at the 
cooling site, as if by Type II supernovae. In this idealized 
immediate feedback heating scenario pH{r, t) is propor- 
tional to q{r)p/tcooi so q is no longer merely a means of 
removing cooled gas, but is directly related to the heat- 
ing. In the 2D version of this heating model, in which q 
is a step function described by Equation (5), we flnd that 
the convcctive turbulence is sufficiently strong to gener- 
ate spontaneous and widely distributed radiative cooling 
within 10 - 20 kpc. The azimuthally averaged gas density 
shows a cyclical variation similar to that described in sec- 
tion 6.2.1. Again, the central temperature is generally too 
high. 

When q is assumed constant (e.g. q = 1), these models 
have some positive features, such as naturally producing 
the observed core in the hot gas density. However, they 
cannot be accepted since multiphase X-ray spectra would 
be expected from the cooling gas and the central masses 
often exceed several 10^ Mq. Type II supernovae have 
never been observed in early type galaxies, so the physical 
origin of the heating process is obscure. 

In summary, by our rigorous standards for success - no 
multiphase gas, agreement with observed profiles, no large 
masses cooling close to the center, etc. - all of the flow 



models we have discussed, as well as many others that we 
have calculated, have failed to match the observations. 

7. MODELS WITH NON- THERMAL PRESSURE 

As we discussed in the Introduction, spatially extended 
non-linear thermal instabilities arise in cooling flows when- 
ever the thermal gas is globally compressed, V • u < 0, or 
9^(r^M)/(9^r < in ID. For a variety of centrally heated 
flows we have shown that runaway radiative cooling is ex- 
pected whenever the normally infiowing gas is decelerated 
and compressed by heated gas closer to center of the flow. 
A similar deceleration can be produced in the absence of 
heating by an additional non-thermal pressure (gradient) 
which is dynamically important within some radius in the 
flow. The non-thermal pressure Pnt could be due to cosmic 
rays or magnetic energy density perhaps associated with 
a central AGN. For simplicity, and in keeping with our 
purely phenomenological approach, we consider two pres- 
surized flows in which the non-thermal pressure is assumed 
to scale with the thermal pressure P: 

Pnt = Pir/rnt)-'. (11) 
The coefficient of P is designed so that dPnt / dr is steeper 
than dP/dr. 

The two columns in Figure 11 illustrate the fiow patterns 
at time t„ in isolated ID galactic flows with r„t = 1 and 2 
kpc respectively. With r„( = 1 kpc the non-thermal pres- 
sure flattens the inner density profile somewhat like the 
NGC 4472 observations. Because the central flow is par- 
tially supported by the additional non-thermal pressure, 
both the thermal pressure and temperature are also lower 
at small r. But when r„t = 1 kpc the non-thermal pres- 
sure gradient is insufficient to keep large masses of gas from 
cooling within 100 pc, M^ooiir < 100 pc,t„) = 1.7 x 10^" 
Mq. However, if the non-thermal pressure is larger, 
r„t = 2 kpc, gas slowly flows out from the core and a 
ID drip forms at r = 0.6 kpc. In this ID flow all the gas 
cools at or near the drip, r = 0.5 — 2.5 kpc, not at the 
origin. Note that Pnt/P must be quite large in this model 
before gas is forced to cool away from the origin. Never- 
theless, neither central or extended cooling is supported 
by XMM observations. 

8. CONCLUSIONS 

For many years it has been assumed that cooling flows 
arc inherently inhomogcneous (multiphase), so that ther- 
mal instabilities develop throughout extended regions of 
the fiow. While there is little or no support for localized 
cooling from linear perturbation studies, by assuming dis- 
tributed cooling dropout it has been possible to consume 
hot gas in a large volume that would otherwise have fiowed 
to the very center and cooled there, generating undesirable 
central mass concentrations and unrealistically peaked X- 
ray emission and thermal temperatures. If this multiphase 
hypothesis is correct, some spectral signature of gas cool- 
ing to low temperatures is expected, such as a broadened 
FeL emission feature. 

It is noteworthy therefore that the first XMM observa- 
tions of elliptical galaxies and clusters have failed to find 
spectral evidence for locally multiphase fiows, even in the 
bright cores of galaxies like M87 (Molendi & Pizzolato 
2001) and NGC 4636 (Xu et al. 2001). These observa- 
tions conflict with FUSE detections of OVI emission from 
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gas at T ^ 3 X 10"'' K observed in several bright elliptical 
galaxies. The OVI emission is a clear indication that gas is 
cooling, although at a somewhat lower rate than expected 
(Bregman, Miller & Irwin 2001). Nevertheless, the null 
XMM observations have revived interest in heating mech- 
anisms that have been proposed to keep gas from cooling 
to low temperatures. If some form of contirnions or inter- 
mittent heating were present, it is claimed, our failure to 
observe cooling and cooled gas with XMM could be easily 
understood. The limited temperature variations that arc 
observed at any projected radius could be explained sim- 
ply by single-phase temperature variations viewed along 
the line of sight. 

To explore this suggestion in more detail we have com- 
puted a large number of cooling flows that arc heated by 
feedback when gas flows into the center. This kind of 
feedback is ideally suited to keep gas from cooling near 
the center of the flow where the gas is densest and radia- 
tive losses are maximal. The heating scenarios we adopt 
are deliberately kept as simple as possible with a mini- 
mum number of free parameters, even though they are 
unphysical in some respects. Our approach is phenomeno- 
logical, with little or no attention to specific (poorly un- 
derstood) microphysical heating processes. We have also 
investigated how thermal conduction and convection can 
share this heating with larger galactic volumes. 

The criteria that we impose for successfully heated mod- 
els arc exacting: (1) the density and temperature profiles, 
n{r) and T{r), at time tn = 13 Gyrs must not differ greatly 
from those observed in typical massive elliptical galaxies 
such as NGC 4472, (2) little or no gas must be currently 
cooling to low temperatures (i.e. <g; 1 keV) anywhere in 
the flow, M ^ 0.1 Mq yr-i for NGC 4472, and (3) the 
total mass of gas that has cooled near the galactic cen- 
ter in the last ~ 10 Gyrs must not exceed the total mass 
typically indicated by observed stellar velocity dispersions. 

After examining a large number of heated cooling flows, 
much larger than we have described here, we conclude that 
it is very difficult to devise any combination of ad hoc 
heating, thermal conduction or convection that produces 
reasonable density and temperature proflles and keeps gas 
from cooling at rates much less than that implied by the 
radiative emission of X-rays, 

which is evaluated for the temperature and bolometric lu- 
minosity of NGC 4472, T « 1.3 x 10^ K and L^^boi « 
7.2 X lO'^^ erg s~^. This result is illustrated most clearly 
in Figure 5, showing that neither the total mass of cooled 
gas nor the present cooling rate M(t„) — dMcooi/dt w 1.2 
Mq yr~^ is appreciably reduced in any of the heated cool- 
ing flows we consider, many of which are extreme in un- 
physical ways. Even continuous heating that exactly bal- 
ances radiative cooling fails spectacularly to agree with 
observations. The central temperatures in heated flows is 
typically too large and the central temperature derivative 
dT/dr is often negative, unlike all well-observed flows. As 
a consequence of these high temperatures, the gas den- 
sities are too low in the center and (often) throughout 
the entire flow. When heating lowers the gas density, we 
find that Type la supernova raise the hot gas iron abun- 
dance far beyond observed values; thus unrealistically high 



iron abundances are an additional serious difficulty with 
heated cooling flows that we have not emphasized here. 
Even when the heating parameters are adjusted for the 
best (but nevertheless unacceptable) fit to the observa- 
tions, the solutions are found to be very sensitive to the 
heating parameters. An undesirable degree of fine-tuning 
appears to be essential for any flow influenced by heating. 

We conclude that the failure of XMM to detect cooling 
gas at the expected rate cannot be easily understood in 
terms of some as yet unspccifled heating mechanism that 
balances or exceeds the X-ray luminosity. 

While the XMM results are inconsistent with multiphase 
flows that cool, optical (Ha -|- [Nil]) and UV line emission 
have been used to support radiative cooling. We estimate 
from the FUSE flux for the OVI doublet in NGC 4636 
(Bregman, Miller & Irwin 2001) that the OVII 21.61(0.574 
keV) line should be about 10 percent as luminous as the 
strong OVIII Lya line detected by Xu et al. (2001) using 
XMM/RGS. Yet Xu et al. were unable to detect the OVII 
line expected from gas cooling near T ~ 2 x 10^ K and 
we have no explanation for this. Nevertheless, our con- 
clusion here remains: if emission from low or intermediate 
temperature ions is absent, this is unlikely to be explained 
by some hypothetical source of heating from the galactic 
center. If the heating is insufficient to drive a global galac- 
tic wind, our heated flow calculations show convincingly 
that the heating actually stimulates gas to cool in a large 
volume of the galaxy. If the heating rate is large enough 
to drive a wind, the gas density (and L^) will be too low 
everywhere, typically by several orders of magnitude. 

On the positive side, we have demonstrated here how gas 
can cool many kpcs from the galactic center when there is 
sonic; cc^iitral heating or non-thermal pressure. Non-linear 
thermal instabilities develop naturally and spontaneously 
in the compressed region of the flow where the normal cool- 
ing inflow from large radii confronts slowed or outflowing 
gas in the heated or pressurized core. Strong local cool- 
ing occurs even when the flow velocities are subsonic. The 
ID representations of this instability, which we refer to as 
"drips" , are unphysical since they are Rayleigh- Taylor and 
convectively unstable. If the cooled gas is removed, spher- 
ically symmetric ID drips usually become quasi-stationary 
at late stages of evolution. If the cooled gas is not removed, 
drips can evolve into massive, high density shells of cooled 
gas that ultimately freefall supersonically into the center of 
the flow, causing unphysical heating and shock reflection. 

It is gratifying therefore that cooling instabilities form 
at nearly the same radius in ID and more realistic 2D 
flows. In 2D flows random, subsonic (convective and tur- 
bulent) velocities induce cooling over a large galactic vol- 
ume. When convective eddies collide, nonlinear compres- 
sions lead to runaway thermal instabilities. The cooling 
is more chaotic in 2D and coherent inward-moving drip 
waves do not form. In principle, gas can also cool at inter- 
mediate radii in unheated galactic cooling flows provided 
the galaxy is rotating (Brighenti & Mathews 1997b, 2000b; 
D'Ercole, &: Ciotti 1998), in this case an extended disk of 
cooling gas is formed, but there is no observational evi- 
dence for this in the X-ray images. We have shown that 
cooling disks can be avoided if the flow is turbulent or 
convective (Brighenti & Mathews 1900) and there is con- 
siderable observational evidence for this (e.g. Caon et al. 
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2000). 

Nonlinear cooling can also be stimulated by an addi- 
tional non-thermal pressure near the center of the cooling 
flow. This pressure could be due to magnetic field stresses, 
cosmic rays or both. Seed magnetic fields ejected from 
mass-losing stars can be amplified by the subsonic turbu- 
lence that results from Type la supernovae or by the ran- 
dom motions observed in cooling flows (Caon et al. 2000). 
The fleld amplification due to this turbulent dynamo satu- 
rates when the magnetic energy density becomes compara- 
ble with that of the subsonic turbulence, which is generally 
less than the thermal energy density (pressure) in the hot 
gas. However, as the field is advected toward the center of 
the cooling flow, B^/8tt grows by adiabatic compression 
and may exceed the gas pressure; at the center (Soker & 
Sarazin 1990; Mathews & Brighenti 1997). Therefore, a 
large non-thermal magnetic pressure could be a natural 
consequence of cooling flows even in the absence of AGN. 
Alternatively, large non-thermal cosmic ray and magnetic 
pressures may be produced by AGN activity associated 
with central black holes. 

Our main conclusions are: 

(1) Cooling flows that are assumed to be centrally heated 
develop nonlinear thermal instabilities on spatial scales 
comparable to that of the hypothetical heating. Thus 
heating from an active nucleus cannot explain the apparent 
absence of multiphase cooling gas in XMM spectra. This 
is particularly true if the heating is large enough to keep 
unrealistic masses of gas from cooling within the central 
kpc. 

(2) Central regions of high nonthermal pressure can also 
stimulate nonlinear thermal instabilities. 

(3) In most or all well-observed cooling flows in galaxy 



groups or clusters the gas temperature has a minimum at 
the center and increases to a maximum typically at several 
stellar effective radii in the group-central E galaxy. The 
higher temperatures in the more distant gas can be sim- 
ply understood in terms of the cosmic accretion of gas onto 
the group dark halo which has a higher virial temperature 
than the central galaxy. As this gas flows inward, its tem- 
perature adapts to the lower virial temperature of the stars 
in the central galaxy, resulting in a positive temperature 
gradient. None of our computed cooling flows with cen- 
tral heating has a monotonically increasing temperature 
profile as observed in this central region. 

(4) If the central heating is strong enough to keep unreal- 

istically large masses of gas from cooling in the central few 
kpc, the global hot gas density profile is usually radically 
altered, and made time dependent. 

(5) Flows heated by central feedback often result in epo- 
sidic cooling that is associated with quasi-cyclic changes 
in the hot gas density profile; at present there is no obser- 
vational evidence for this. 
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Fig. 1. — Gas temperature, density and radial velocity at time tn = 13 Gyrs for ID isolated (solid lines) and accreting (daslied lines) 
galactic cooling flow models with two cooling dropout assumptions. Left column: Flows with no distributed cooling dropout, q = 0. Right 
column: Flows with uniformly distributed cooling, q = 1. The temperatures are shown as functions of physical radius, not averaged in 
projection along the line of sight. Density and temperature observations for NGC 4472 are Einstein data (Trinchieri, Fabbiano, & Canizajres 
1986) (filled circles) and ROSAT data (Irwin & Sarazin 1996) (open circles). 

Fig. 2. — Gas temperature, density, specific entropy and velocity at time t„ = 13 Gyrs for four ID heated galactic cooling flows. Each 
panel shows flows for isolated (solid lines) and accreting (dashed lines) galaxies. The temperatures are shown as functions of physical radius, 
not averaged in projection along the line of sight. Left Column: A galactic cooling flow model in which cooling is locally balanced by heating. 
Right Column: Comparison of flows heated as described in Eq. (8) with Gh = 1 erg s ~^ and Vfi = 5 kpc. A drip is seen at r s» 7 — 8 
kpc for both flows. Data is the same as in Figure 1. 

Fig. 3. — The integrated mass of cooled gas Mcool(''') at various times for the heated cooling flow in the isolated E galaxy shown in the 
right column of Figure 2. Note the M^ggi(r) plot is lineaj:, not logarithmic. 

Fig. 4. — Four panels show the 2D flow with Gaussian heating (G^ = 1, = 5 kpc) at time tn = 13 Gyrs. Upper Left: Gas velocity and 
density. Contours show log p(R, z). The velocity arrow at the upper right represents 200 km s~^. Upper Right: Density contours of cooled 
gas. Lower Left: Azimuthally averaged hot gas density shown with the same data as in Figure 1. Lower Right: Azimuthally averaged cooled 
gas density shown at three times. 

Fig. 5. — Plots of the growth of cooled gas mass M^aoi{t) for a large number of heated cooling flows. All flows arc for the "isolated" galaxy. 
Flows with Gaussian heating are shown with the values of the two heating parameters in Equation (8), Gf^ and r^. Gaussian- heated flows 
computed in 2D are labeled; all other flows are in ID including the two flows with conduction and convection. Flows heated as described by 
Equation (9) are labeled CH; weak and strong heating refer respectively to C^Tc = 1.6 x 10*^ and 1.6 x 10** erg cm-^ g^^ s^^. Top Panel: 
Mcool(t) for four Gaussian heated flows all with Gh = 1 erg gm~^ s~^ and r/j = 5 kpc) compared with a flow with no heating. Middle 
Panel: M^„„i(t) for three Gaussian heated flows and two CH flows. For the (G/j,r/j)= (0.5,20) flow the cooled mass within 10 kpc is shown 
separately. Bottom Panel: M^ggi(t) for three 2D heated flows. 

Fig. 6. — The density cycle during a cooling event. Azimuthally averaged hot gas density profiles for the 2D flow heated with (G^, r/^) = 
(1, 5 kpc) are shown at several times preceding, during and following a cooling episode. Times in Gyrs label each density profile. The times 
and total masses of cooled gas are: (t Gyrs; M^gai/lO^" Mq) = (8.0; 3.185), (9.1; 3.19), (9.4; 3.194), (9.7; 3.283), (9.8; 3.618), (10; 3.660) and 
(11; 3.685). Data is the same as in Figure 1. 

Fic;. 7. — Gas temperature, density, specific entropy and velocity at time t„ = 13 Gyrs for four Gaussian heated cooling flow models. 
Flows for isolated and accreting galaxies are shown with solid lines and dashed lines respectively. Gas temperatures are shown as functions 
of physical radius, not averaged in projection along the line of sight. First Column: Flow for (Gf^, rji) = (0.5, 5 kpc). Second Column: Flow 
for (Gh, f fe) = (2, 5 kpc). Third Column: Flow for (G^, r^) = (0.5, 20 kpc). Fourth Column: Flow for (G^, r^) = (20, 1 kpc). Data is the 
same as in Figure 1. 

Fig. 8. — Gas temperature, density, specific entropy and velocity at time tn = 13 Gyrs for two cooling flow models with Compton-like 
heating described by Equation (9). Flows for isolated and accreting galaxies are shown with solid lines and dashed lines respectively. Gas 
temperatures are shown as functions of physical radius, not averaged in projection along the line of sight. Left Column: Weakly heated flow 
with C^Tq = 1.6 X 10*^ erg cm^ g"^ s~^. Right Column: Strongly heated flow with GfiTc = 1.6 x 10**. Data is the same as in Figure 1. 

Fig. 9. — Representative flows that are continuously heated by hypothetical stellar sources, all shown at time tn = 13 Gyrs. This type 
of SNH heating is implemented by artificially increasing the Type la supernova rate SNu(tn), each supernova releasing 10®^ ergs. Values of 
SNu(t„) in supernovae during 100 years per 10^" label each curve in the top panel and apply for lower panels in each column. Unlike 
previous figures, isolated and accretion flows are shown in separate columns. Data is the same as in Figure 1. 

Fig. 10. — Three 2D flows at time t„ = 13 Gyrs for an isolated galaxy. The uppermost panels show the gas density and velocity; the 
velocity arrow at the upper right represents 200 km s~^. The second row of panels shows density contours of cooled gas. The third and 
bottom rows of panels show the azimuthally averaged hot gas and cooled gas densities respectively. Left Column: Gaussian heated flow with 
(G/j, r/j) = (20, 1 kpc). Central Column: Flow weaJjly heated with Equation (9) using C^Tc = 1.6 x 10*^ erg cm^ g~^ s~^. Right Column: 
Strongly heated flow with GhTc = 1.6 x 10**. Data is the same as in Figure 1. 

Fig. 11. — Gas temperature, density, pressure and velocity at tn = 13 Gyrs for two unheated isolated ID flows with additional non-thermal 
pressure described by Equation (11). Left Column: Flow for rnt = 1 kpc. Right Column: Flow for rnt = 2 kpc. The third pair of panels from 
the top compare the profiles of thermal pressure (solid lines) and non-thermal pressure (dashed lines). 
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APPENDIX 
APPENDIX 

We describe how the ID Gaussian heated flows with G/i = 1 erg gm"^ s~^ and Vh = 5 kpc are influenced by thermal 
conduction and convection. 

Influence of Thermal Conduction 

The conductive heat flux is 

F ^ 

cona — n 

or 

where k = /1. 84 x 10~^T^/^/lnA is the thermal conductivity, which can be reduced by a factor / < 1 to approximate 
the inhibiting eS'ect of tangled magnetic fields. When the conduction becomes saturated, the fiux is replaced with 

Fcond,sat = QAfnkT{2kT/TTmef'^. 

Because k is very small at low temperatures, thermal conduction can usiially be ignored for cooling flows in galaxies and 
galaxy groups where T < 1 keV. However, in the flow near and inside drips the temperature gradients are very large, 
increasing the likelihood that F^ond 

may alter the solution and possibly the mass dropout stimulated by heating at the 

galactic center. 

The left column of Figure Al shows the flow variables at tn = 13 Gyr for isolated and accreting flows with the 
conductivity set at the full Spitzcr value, / = 1. The drip that formerly appeared at 7 kpc (right column of Figure 2) has 
(temporarily) disappeared at time tn when there are no thermally unstable compression regions. While the flows with 
conduction are smoother overall, there is again an undesirable negative temperature gradient in the central region. 

A closer look at these thermally conducting gas flows reveals that many drips occurred at t < tn and a large amount of 
gas cooled (dropped out) at large radii. A large mass of gas, ^ 2.0 x 10^° Mq, also cools in r < 100 pc, where we assume 
g = 1 to absorb the cooled gas. Cooling dropout that occurred before t — 3 Gyrs, when Mcooi ~ 10^° Mq, may have 
formed into a second generation of old stars. By time tn, Mcooi = 3.7 x 10^° (7.9 x 10^°) Mq has cooled in the isolated 
(accretion) flow. Most of the cooled gas is concentrated near r « 7 kpc where the drip appeared without conduction. 

At late times, episodic drips occur in these flows separated by '--^ 2 Gyr periods of relative quiescence, resulting in a 
quasi-cyclic variation in the global gas density as discussed in §6.2.1. Just following a period of intense spatially distributed 
mass dropout in drip waves, the hot gas density is lower and the drip cooling shuts down temporarily. During the quiescent 
period that follows, the gas density slowly increases due to stellar mass loss until it becomes high enough to stimulate 
another episode of drip cooling. The high gas density in r ;^ 10 kpc visible in the accretion flow with conduction occurs 
at a moment jiist prior to a period of more intense cooling. Two flows that appear to be quite different at time i„, such 
as the isolated and accretion flows with conduction in Figure Al differ largely because they are momentarily at different 
phases of the cooling cycle. 

Since neither of the flows with thermal conduction shown in Figure Al is currently undergoing drip cooling, they would 
show no spectral evidence for multiphase gas at the present time. Nevertheless, at any particular time about half of all 
E galaxies with this type of flow should be actively cooling and exhibit multiphase spectra. In addition, the negative 
temperature gradients would also disqualify both solutions with full Spitzer (/ = 1) conduction. When these calculations 
were repeated with a reduced conductivity, / = 0.1, which may be more astrophysically plausible, the drip structure and 
cooling dropout become almost identical to that of the / = case shown in the right column of panels in Figure 2. We 
conclude that thermal conduction is unable to significantly reduce the distributed cooling and total mass of cooled gas in 
flows heated by the Gaussian H{r,t). 

Influence of Convection 

Since the entropy gradient is negative near drip waves, convection is expected. In ID flow the convective flux is 

'dlogT\ 7-1" 

where cp = j/k/[{-f — l)^mp], g is the local acceleration of gravity, and Hp = P/pg is the pressure scale height 
(Kippenhahn & Weigert 1990). This stellar expression for Fconv is appropriate since cooling flows are very nearly in 
hydrostatic equilibrium. We have explored various values of the mixing length £ with similar results. In the right column 
of Figure Al we illustrate the same Gaussian heated flow at time t„ as in the right column of Figure 2 but with convective 
heat flux (and Fcond = 0) using £ = max(2r/3, ZHp/A), chosen to maximize the effect of convection. The results are very 
similar to the conductive case. Frequent episodic drips still occur, producing a maximum in the density of cooled gas at 
r ~ 5 kpc at i„. The total mass of cooled gas does not change appreciably when convection is included. 

In the right columns of Figures 2 and Al we also show the entropy proflles S{r) = {cp/^)\n.{T /n^/^) at time tn for 
Gaussian heated flows with and without convection {Gh = 1 erg gm~^ s~^ and r/t = 5 kpc). As expected, the effect of 
convection is to flatten the steep entropy profile in the drip structure. However, for the mixing length parameters we 
consider the entropy gradient is still mildly negative (in 4 < r < 20 kpc) even in the presence of convection. 

The convective flows shown in Figures Al - with its unrealistically peaked temperature profile and expected multiphase 
spectrum - underscore our conclusion: spherically symmetric heated cooling flows cannot be brought into agreement with 
XMM observations by including convection. 
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Fig. A12. — Gas temperature, density, specific entropy and velocity at time tn = 13 Gyrs for ID Gaussian heated galactic cooling flows 
with Gh = 1 erg g~^ s ~^ and r/j = 5 kpc. Each panel shows flows for isolated (solid lines) and accreting {dashed lines) galaxies. The 
temperatures are shown as functions of physical radius, not averaged in projection along the line of sight. Left Column: Gaussian heated 
flows identical to those in the second column of Figure 2 but including thermal conduction at the full Spitzer value (/ = 1). Right Column: 
Flows identical to those in the second column of Figure 2 but including convection. Data is the same as in Figure 1. 
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